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ABSTRACT: 

In this note we initiate the probabilistic study of the critical points of polynomials of 
large degree with a given distribution of roots. Namely, let / be a polynomial of degree 
n whose zeros are chosen IID from a probability measure /x on C. We conjecture that 
the zero set of /' always converges in distribution to /i as n — )■ oo. We prove this for 
measures with finite one-dimensional energy. When /i is uniform on the unit circle 
this condition fails. In this special case the zero set of /' converges in distribution 
to that the IID Gaussian random power series, a well known determinantal point 
process. 
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1 Introduction 



Since Gauss, there has been considerable interest in the location of the critical points 
(zeros of the derivative) of polynomials whose zeros were known - Gauss noted that 
these critical points were points of equilibrium of the electrical field whose charges 
were placed at the zeros of the polynomial, and this immediately leads to the proof 
of the well-known Gauss-Lucas Theorem, which states that the critical points of a 
polynomial / lie in the convex hull of the zeros of / (see, e.g. |Mar49t Theorem 6,1]). 
There are too many refinements of this result to state. A partial list (of which several 
have precisely the same title!) is as follows: |Mar83t IDro89[ IStoQGj IPaw98t ITarSG 



[Rih72l IGRR691 [Joy69| ldBS47l IdBlGl IBra31[ IMahGll [Dm98l IMalOSl ICM04 
ISenlOt ISenOlj ). Among these, we mention two extensions that are easy to state. 



• Jensen's theorem: if p{z) has real coefficients, then the non-real critical points of 
p lie in the union of the "Jensen Disks" , which are disks one of whose diameters 
is the segment joining a pair of conjugate (non-real) roots of p. 

• Marden's theorem: Suppose the zeroes zi,Z2, and of a third-degree polyno- 
mial p{z) are non-coUinear. There is a unique ellipse inscribed in the triangle 
with vertices zi, Z2, Z3 and tangent to the sides at their midpoints: the Steiner 
inellipse. The foci of that ellipse are the zeroes of the derivative p'{z). 



There has not been any probabilistic study of critical points (despite the obvious 
statistical physics connection) from this viewpoint. There has been a very extensive 
study of random polynomials (some of it quoted further down in this paper), but 
generally this has meant some distribution on the coefficients of the polynomial, and 
not its roots. Let us now define our problem: 

Let /i be a probability measure on the complex numbers. Let {X„ : n > 0} 
be random variables on a probability space {Q, J-", P) that are IID with common 
distribution /i. Let 

n 

fn{z) ■.= 1[{Z-X,) 
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be the random polynomial whose roots are Xi, . . . ,X„. For any polynomial / we 
let 2{f) denote the empirical distribution of the roots of /, for example, Z{fn) = 

The question we address in this paper is: 
Question 1.1. When are the zeros of f'^ stochastically similar to the zeros of fn? 

Some examples show why we expect this. 

Example 1.1. Suppose /i concentrates on real numbers. Then /„ has all real zeros 
and the zeros of interlace the zeros of /„. It is immediate from this that the 
empirical distribution of the zeros of converges to /i as n — )■ oo. The same is true 
when /i is concentrated on any affine line in the complex plane: interlacing holds and 
implies convergence of the zeros of to /i. Once the support of /i is not contained in 
an affine subspace, however, the best we can say geometrically about the roots of 
is that they are contained in the convex hull of the roots of /„; this is the Gauss-Lucas 
Theorem. 

Example 1.2. Suppose the measure /i is atomic. If /i(a) = p > then the multiplicity 
of a as a zero of /„ is n{p + o{l)) . The mulitplicity of a as a zero of f^ is one less 
than the multplicity as a zero of /„, hence also n{p + o(l)). This is true for each of 
the countably many atoms, whence it follows again that the empirical distribution of 
the zeros of f^ converges to /i. 

Atomic measures are weakly dense in the space of all measures. Sufficient conti- 
nuity of the roots of /' with respect to the roots of / would therefore imply that the 
zeros of always converge in distribution to /i as n — )■ oo. In fact we conjecture this 
to be true. 

Example 1.3. Our first experimental example has the roots of / uniformly dis- 
tributed in the unit disk. In the figure, we sample 300 points from the uniform 
distribution in the disk, and plot the critical points (see Figured]). The reader may 
or may not be convinced that the critical points are uniformly distributed. 
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Figure 1: Critical points of a polynomial whose roots are uniformly sampled inside 
the unit disk. 




Figure 2: Critical points of polynomial whose roots are uniformly sampled on the 
unit circle. 
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Example 1.4. Our second example takes polynomials with roots uniformly dis- 
tributed on the unit circle, and computes the critical points. In Figure [2] we do 
this with a sample of size 300. One sees that the convergence is rather quick. 

Remark 1.1. The figures were produced with Mathematica. However, the reader 
wishing to try this at home should remember to do all the computations to very high 
precision, otherwise Mathematica (and Matlab and R) return garbage answers, since 
they all use a rather primitive method of computing zeros of polynomials. 

Conjecture 1.1. For any fi, as n oo, Z{f') converges weakly to fi. 

There may indeed be such a continuity argument, though the following counterex- 
ample shows that one would at least need to rule out some exceptional sets of low 
probability. Suppose that f{z) = — 1. As n — )• oo, the distribution of the roots 
of / converge weakly to the uniform distribution on the unit circle. The roots of 
however are all concentrated at the origin. If one moves one of the n roots of /„ along 
the unit circle, until it meets the next root, a distance of order 1/n, then one root of 

zooms from the origin out to the unit circle. This shows that small perturbations 
in the roots of / can lead to large perturbations in the roots of /'. It seems possible, 
though, that this is only true for a "small" set of "bad" functions /. 

1.1 A little history 

This circle of questions was first raised in discussions between one of us (IR) and 
the late Oded Schramm, when IR was visiting at Microsoft Research for the auspi- 
cious week of 9/11/2001. Schramm and IR had some ideas on how to approach the 
questions, but were somewhat stuck. There was always an intent to return to these 
questions, but Schramm's passing in September 2008 threw the plans into chaos. We 
(RP and IR) hope we can do justice to Oded's memory. 

These questions are reminscent of questions of the kind often raised by Herb Wilf, 
that sound simple but are not. This work was first presented at a conference in Herb's 
honor and we hope it serves as a fitting tribute to Herb as well. 
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2 Results and Notations 



Our goal in this paper is to prove cases of Conjecture 11.11 



Definition 2.1. We definite the p- energy of fi to he 




1 



i/p 



z — w\P 



Since in the sequel we will only be using the 1-energy, we will write S for Si. 
By Fubini's Theorem, when has finite 1-energy, the function defined by 



Then 2{f^) converges in distribution to ^ as n —t- oo. 

A natural set of examples of fi with finite 1-energy is provided by the following 
observation: 

Observation 2.3. Suppose Q G C has Hausdorff dimension greater than one, and /i 
is in the measure class of the Hausdorff measure on Vl. Then fi has finite 1-energy. 

Proof: This is essentially the content of |Fal03] [Theorem 4.13(b)]. □ 

In particular, if fi is uniform in an open subset (with compact closure) of C, its 
1-energy is finite. 

A natural special case to which Theorem 12.21 does not apply is when /i is uniform 
on the unit circle; here the 1-energy is just barely infinite. 




is well defined and in L^{fi). 



Theorem 2.2. Suppose jj, has finite 1-energy and that 



/i{^:V;(z) = 0} = 0. 



(2.1) 
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Theorem 2.4. // fi is uniform on the unit circle then Z{fn) converges to the unit 
circle in probability. 

This resuh is somewhat weak because we do not prove Z{fn) has a hmit in dis- 
tribution, only that aU subsequential hmits are supported on the unit circle. By the 
Gauss-Lucas Theorem, all roots of /„ have modulus less than 1, so the convergence 
to /i is from the inside. Weak convergence to /x implies that only o{n) points can be 
at distance 0(1) inside the cirle; the number of such points turns out to be 0(1). 
Indeed quite a bit can be said about the small outliers. For < p < 1, define 
Bp := {z : \z\ < p}. The following result, which implies Theorem 12. 4[ is based on a 
very pretty result of Peres and Virag |PV05t Theorems 1 and 2] which we will quote 
in due course. 

Theorem 2.5. For any p G (0, 1), as n ^ oo, the set Z{gn) H Bp of zeros of Qn on 
Bp converges in distribution to a determinantal point process on Bp with the so-called 
Bergmann kernel 7r^-'^(l — ZiZjY . The number N{p) of zeros is distributed as the sum 
of independent Bernoullis with means p^^ , 1 < k < oo. 



2.1 Distance functions on the space of probability measures 

If p and V are probability measures on a separable metric space then the Pro- 
horov^l distance \p — z/|p is defined to be the least e such that for every set A, 
p{A) < z^(v4'^) + e and ^{A) < piA"^) + e. Here, A"^ is the set of all points within 
distance e of some point of A. The Prohorov metric metrizes convergence in distribu- 
tion. We view collections of points in C (e.g., the zeros of /„) as probability measures 
on C, therefore the Prohorov metric serves to metrize convergence of zero sets. The 
space of probability measures on S, denoted V{S), is itself a separable metric space, 
therefore one can define the Prohorov metric on V{S), and this metrizes convergence 
of laws of random zero sets. 

The Ky Fan metric on random variables on a fixed probability space will be of 
some use as well. Defined by K{X,Y) = inf{e : F{d{X,Y) > e) < e}, this metrizes 



'Also known as the Prokhorov and the Levy-Pro(k)horov distance 
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convergence in probability. The two metrics are related (this is Strassen's Theorem): 

l/i - v\p = mi{K{X, F) : X ~ /i, F ~ z/} . (2.2) 

A good reference for the facts mentioned above is [Hofoej^ We will make use of 
Rouche's Theorem. There are a number of formulations, of which the most elementary 
is probably the following statement proved as Theorem 10.10 in |BN82j . 

Theorem 2.6 (Rouche). /// and g are analytic on a topological disk, B, and \g\ < |/| 

on OB , then f and f + g have the same number of zeros on B. □ 



3 Proof of Theorem [23 



We begin by stating some lemmas. The first is nearly a triviality. 
Lemma 3.1. Suppose fi has finite 1-energy. Then 



t-p(^|Xo-Xi|<i) ^0. 

(a) for any C > 0, 

P I min |X,- - Xn+i\ < — 1 ^ 

\l<3<n n ' 

Proof: For part (i) observe that limsupt ■ P(|Xo - Xi\ < 1/t) < 21imsup2-'' 
P (|Xo — Xi\ < 2^^) as if: goes over reals and j goes over integers. We then have 



oo > 

= El/|Xo-Xi| 

> ^E5^2n|x„_x,,<2- 



2 



IJ2^^¥{\Xo-X,\<2-^) 



2 

3 



^available on line at ,http : / / epub . oeaw . ac . at/0xclaa500d_ 0x00239061 . pdf 
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and from the finiteness of the last sum it follows that the summand goes to zero. 
Part (m) follows from part (i) upon observing, by symmetry, that 

P ( min |X,- - Xn+i\ < —] < ( IXq - XJ < - ) . 
\i<j<n n J \ n J 

□ 

Define the n*^ empirical potential function „ by 



1 " 1 

n ^—^ r — Y ■ 



n z — Xj 



which is also the integral in w of l/{z — w) against the measure Z{fn). Our next 
lemma bounds V^„(-2) on the disk B := Bc/n{^n+i)- 



Lemma 3.2. For all e > 0, 



¥(snp\V;^Jz)\>en) ^0 
\zeB J 



as n oo. 



Proof: Let G„ denote the event that mini<j<„ |Xj — X^+il > 2C/n. Let Sn '■- 
sup^eB l^,n(^)l- We will show that 



E^„1g„ = o{n) 



(3.1) 



n — > oo. By Markov's inequality, this implies that P(S'„1g„ > ^n) — )■ for all e > as 
n — )■ oo. By part (ii) of Lemma [3.11 we know that P(G'„.) ~^ 1, which then establishes 
that W{Sn > €n) — 0, proving the lemma. 



In order to show (13.11) we begin with 

n 

-E 



< 



1^ 1 

\z - X, 



12 ■ 



Therefore, 



-E 



+1 



C/n) 



< 



n 

-E 



X„+i — Xj 
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and by symmetry, 

ESnlG„ < 4E|Xo - X,\-H\x,-X,\>2C/n < 4E — — A — . 

V 1-^0 ~ ^j\ / 

We may then integrate by parts, obtaining 

E5„1g„ < 8tP (^^— ^ >tj dt. 

The integrand goes to zero as n — t- oo by part (z) of Lemma 13.11 It follows that the 
integral is o{n), proving the lemma. □ 

Define the lower modulus of V to distance C/n by 

inf \V,,nH\ . 

w.\w—z\<C /n 

This depends on the argument /i as well as C and n but we omit this from the 
notation. 

Lemma 3.3. Assume fi has finite 1-energy. Then as n oo, the random variable 
V_^{Xn+i) converges in probability, and hence in distribution, to \Vf^{Xn+i)\- 

In the sequel we will need the Glivenko-Cantelli Theorem |Dur04t Theorem 1.7.4]. 
Let Xi . . . ,Xn, ... be independent, identitically distributed random variables in M 
with common cumulative distribution function F. The empirical distribution function 
Fn for Xi, ... , Xn is defined by 



1 " 

Fn{x) = - ^I(-oo,x]{Xi 



n . 
1=1 



where Ic is the indicator function of the set C. For every fixed x, is a sequence 

of random variables, which converges to F{x) almost surely by the strong law of large 
numbers. Glivenko-Cantelli Theorem strengthen this by proving uniform convergence 
of Fn to F. 



Theorem 3.4 (Glivenko-Cantelli). 



l-^n — -^lloo = sup \Fn{x) — F{x)\ — > almost surely. 
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The following Corollary is immediate: 
Corollary 3.5. Let f be a bounded continuous function on M. Then 

lim / fdFn = / fdF, almost surely. 



Another immediate Corollary is: 

Corollary 3.6. With notation as in the statement of Theorem 3.4, the Prohorov 
distance between F^ and F converges to zero almost surely. 

Proof of Lemma I3.3t It is equivalent to show that - |V"^(X„+i)| in 
probability, for which it sufficient to show 

sup \V^,n{u) - V^{Xn+i)\ ^ (3.2) 

u&B 

in probability. This will be shown by proving the following two statements: 

sup \V^,n{u) - V^^n{Xn+i)\ in probability ; (3.3) 

|V;.,n(X„+i) - V;.(X„+i)| ^ O in probability . (3.4) 

The left-hand side of (13.31) is bounded above by (C/n) sup„g^ | V^^^^ {u)\. By Lemma [3^ 
for any e > 0, the probability of this exceeding Ce goes to zero as n — )■ oo. This es- 
tablishes (13. 3p . 

For (13. 4p we observe, using Dominated Convergence, that under the finite 1-energy 
condition, 

S^{fi):= / / -l\^_^\-i>K dij,{z) djj,{w) ^ 



\z — w\ 



as — )■ oo. Define 0^'^ by 



\z — w\ 



z — w max{ I z — w 1 , 1/ K} 

in other words, it agrees with l/{z — w) except that we multiply by a nonegative real 
so as to truncate the magnitude at K. We observe for later use that 

1 

\z — w\ \z — w\ 



0^'^(«;) 



< ■: rl|2-tu|-i>x 
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so that 




[w - 



\z — w\ 



dn{z) dfi{w) < S"^ ifi) ^ . (3.5) 



We now introduce the truncated potential and truncated empirical potential with 



respect to 0^'' 



We claim that 



Indeed, 



E |V;^(X„+i) - V;(X„+i)| < S^ifi) . (3.6) 



V,{X^^,) - V;^(X„+i) = j (73^ - ^''''(^n+i)) dii{z) 

so taking an absolute value inside the integral, then integrating against the law of 
Xn+i and using fl3.5p proves (13.61) . The empirical distribution „ has mean and is 
independent of therefore the same argument proves 

E |v;^„(x„+i) - v;,„(x„,+i)| < ^""(z^) (3.7) 

independent of the value of n. 

We now have two thirds of what we need for the triangle inequality. That is, to 
show (13. 4p we will show that the following three expressions may all be made smaller 
than e with probability 1 — e. 

Vfj,^n{Xn+l) - V^^(X„+i) 

Choosing K large enough so that S^{fi) < e^, this follws for the third of these follows 
by (13.61) and for the first of these by (13.71) . Fixing this value of K, we turn to the 
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middle expression. The function 0^'^ is bounded and continuous. By the Corollary 
l3.5l to the Glivenko-Cantelli Theorem 13. 4[ the empirical law Z{fn) converges weakly 
to fi, meaning that the integral of the any bounded continuous function (p against 
Z{fn) converges in probability to the integral of (p against fi. Setting (p := 0^'^ and 
z := Xn+i proves that \/^(X„+i) — V^(X„_i_i) goes to zero in probability, establishing 
the middle statement (it is in fact true conditionally on Xn+i) and concluding the 
proof. □ 

Proof of Theorem [221 Suppose that V_^iXn+i) > l/C. Then for all w with 

\w — Xn+il < C/n, we have 

j=i J 



and hence 

7? 

I/;HI = n \V,,^{w)\ > nV^{X^+,) > - . 

To apply Rouche's Theorem to the functions 1//^ and z — X^+i on the disk B : = 
Bc/n{.Xn+i) we note that |1//^| < C/n = \z — on dB and hence that the sum 

has precisely one zero in B, call it a„+i. Taking reciprocals we see that a„+i is also the 
unique value 'm. z & B for which fn{z) = —l/{z — Xn+i)- But fn{z) + l/(z — Xn+i) = 
fn+ii^)y whence f^_^_i has the unique zero a„+i on B. 

Now fix any 6 > 0. Using the hypothesis that ^{z : V^{z) = 0} = 0, we pick a 
C > such that P(|V^(X„+i)| < 2/C) < 6/2. By Lemma [331 there is an uq such 
that for all n > uq, 

It follows that the probability that /^^^ has a unique zero a^+i in B is at least 1 — 5 
for n > uq. By symmetry, we see that for each j, the probability is also at least 
1 — 6 that f^_^_^ has a unique zero, call it aj, in the ball of radius C/n centered at Xj; 
equivalently, the expected number of j < n + 1 for which there is not a unique zero 
of fn+i in Bc/n{Xj) is at most 6n for n > no. 

Define Xj to equal aj if f^^-^^ has a unique root in Bc/n{Xj) and the minimum 
distance from Xj to any Xi with i < n+1 and i ^ j is at least 2C/n. By convention, 
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we define Xj to be the symbol A if eitlier of these conditions fails. The values Xj 
other than A are distinct roots of f^^^ and each such value is within distance C/n of 
a different root of fn+i- Using part (ii) of Lemma [STT] we see that the expected number 
of j for which xj = A is o{n). It follows that P(|Z(/„+i) - ZUl^^)\p > 26) ^ as 
n — 7- oo. But also the Prohorov distance between Z{fn+i) and fi converges to zero by 
Corollary 13.61 The Prohorov distance metrizes convergence in distribution and 6 > 
was arbitrary, so the theorem is proved. □ 



4 Proof of remaining theorems 

Proof of Theorem[2]51 Let Q := ^JLqYjz^ denote the standard complex Gaussian 
power series where {Yj{u)} are IID standard complex normals. The results we require 
from |PV05j are as follows. 

Proposition 4.1 ( |P V05] ) . The set of zeros of Q in the unit disk is a determinantal 
point process with joint intensities 

p{zi, ...,Zn) = 7r""det 

The number N{p) of zeros of Q on Bp is distributed as the sum of independent 
Bernoullis with means p^^ , 1 < k < oo. □ 

To use these results we broaden them to random series whose coefficients are 
nearly IID Gaussian. 

Lemma 4.2. Let {gn '■= X]-^o '^"^'^''} ^ sequence of power series. Suppose 

(i) for each k, the k-tuple (a„,i, . . . , ctn.fc) converges weakly as n ^ oo to a k-tuple 
of IID standard complex normals; 

(ii) E|anr| < 1 for all n and r. 
Then on each disk Bp, the set Z{gi) fl Bp converges weakly to Z{Q) fl p. 



1 

(1 - Z.[ZjY 
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Proof: Throughout the proof we fix p e (0, 1) and denote B :— Bp. Suppose an 
analytic function h has no zeros on dB. Denote by | /i| |b the sup norm on functions 
restricted to B. Note that ii hn ^ h uniformly on B then Z{hn) r\ B ^ Z{h) f] B in 
the weak topology on probability measures on B, provided that h has no zero on dB. 
We apply this with h — Q :— Y1^=q^3^^ where {Yj{u)} are IID standard complex 
normals. For almost every cu, h{u:) has no zeros on dB. Hence given e > there 
is almost surely a 5{uj) > such that H^? — < S imphes \2{g) — Z{Q)\p < e. 
Pick So{e) small enough so that F{6{u) < Sq) < e/3; thus Hf? — < Sq implies 
\Z{g) — Z{Q)\ < e for all Q outside a set of measure at most e/3. 



By hypothesis (ii), 



E 



E 

r=k+l 



< 



I-P 



Thus, given e > 0, once k is large enough so that p /(I — p) < e5o(e)/6, we see that 



P 



Or, 

r=k+l 




For such a k{e) also | Xl^fc+i ^-^''1 — By hypothesis (i), given e > and the 

corresponding 5{e) and A;(e), we may choose no such that n > tiq implies that the 
law of (a„i, . . . ,a„fc) is within min{e/3, 5o(e)/(2A;)} of the product of k IID standard 
complex normals in the Prohorov metric. By the equivalence of the Prohorov metric 
to the minimal Ky Fan metric, there is a pair of random variables g and h such that 
g ^ g-n and h Q and, except on a set of of measure e/3, each of the first k coefficients 
of g is within 60/ {2k) of the corresponding coefficient of Q. By the choice of k{e), we 
then have 

n\\~9 

By the choice of Sq, this implies that 



r,, c- X 2e 
h\\B>So)<j. 



F{\Z{~g)-ZCh)\p >€)<€. 

Because g ^ gn and /i ~ ^, we see that the law of Z{gn) n B and the law of Z{Q) n B 
are within e in the Prohorov metric on laws on measures. Because e > was arbitrary, 
we see that the law of Z{gn) n B converges to the law of Z{Q) (IB. □ 
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Proof of Theorem 12.51 Let p < 1 be fixed for tlie duration of this argument and 



denote B := Bp. Let 



9n[Z) 



Because \Xj\ = 1, the rational function l/(z — Xj) = —X^^/{1 — Xj^z) is analytic 
on the open unit disk and represented there by the power series — YlT=o-^i^~^'^'' ■ 
follows that — (?„/ ^Jn is analytic on the open unit disk and represented there by the 
power series —gn{z)/y/n = YlT=o'^nrZ^ where 

a^r = n-'''Y.Xj^-\ 
i=i 

The function —Qn/y/n has the same zeros on B as does the normalization by 
— 1/y/n being inserted as a convenience for what is about to come. 

We will apply Lemma 14.21 to the sequence {(?„}. The coefficients a„j are nor- 
malized power sums of the variables {Xj}. For each r > and each j, the vari- 
able X~^~^ is uniformly distributed on the unit circle. It follows that Eonr = 

and that Ea^^a^ = n~^^^-X^^~^Xj ^ = n~^J2ij^ij = 1- particular, 
Ejcinrl < (E|a„rP)^''^ = 1, satisfying the second hypothesis of Lemma 14.21 For 
the first hypothesis, fix k, let 9j = Arg(Xj), and let v'^-'^ denote the (2A;)-vector 
(cos(^j), — sin(^j), cos(2^^j), — sin(2^^j), . . . , cos{k9j), — sin(A;6'j)); in other words, v*^-^^ 
is the complex fc- vector (X^~^,X~^, . . . ^X''^) viewed as a real (2A;)-vector. For each 

1 < s,t < 2k we have Evi"'''vp'' = {l/2)6ij. Also the vectors {v'^-'^} are independent 
as j varies. It follows from the multivariate central limit theorem (see, e.g., |Dur04t 
Theorem 2.9.6]) that u*^"^ := n"-*^/^ ^"^-j^ v^-'^ converges to 1/a/2 times a standard 
(2fc)-variate normal. For 1 < r < k, the coefficient a„r is equal to Ug^Li + iu'^\ Thus 
{o-nr : 1 < T < /c} couvcrgcs in distribution as n — )■ oo to a A;-tuple of IID standard 
complex normals. The hypotheses of Lemma 14.21 being verified, the theorem now 
follows from Proposition 14.11 □ 
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